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ABSTRACT 

A Roche model for describing uniformly rotating rings is presented and the results are com- 
pared with numerical solutions to the full problem for poly tropic rings. In the thin ring limit, 
the surfaces of constant pressure including the surface of the ring itself are given in analytic 
terms, even in the mass-shedding case. 
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1 INTRODUCTION 

A uniformly rotating star with a sufficiently soft equation of state 
can b e described approximately using the Roche model iRochd 
In this model the matter is treated as a test fluid in a 
1/r-potential, i.e. one considers the whole mass of the star to 
be concentrated at the centre. By doing so, self-gravitating ef- 
fects of the outer mass shells are neglected. Such Roche mod- 
els have_been_considered, espec ia lly in the m ass-shedding limit, 
by IZel'dovicF & Noviko\l h97lh : IShapiro & S hibata (2003). For 
polytropes with indices ti > 2.5 this approximation yields results 
that differ by less t han about a percent from their correct values 
jMeinel et al.ll2008[) . Whereas most analytical solutions to figures 
of equilibrium describe bodies with constant mass density (e.g. the 
Maclaurin spheroids, lacobi ellipsoids, etc.), the Roche model is 
very useful because it is applicable to non-homogeneous matter. 

Aside from spheroidal figures of equilibrium, there exist 
configurations with toro idal topology . Suc h rings have been 
studi ed both analytically (Kowalewskv 1885; Poincare 1885a, b e; 
iDvson 1892. 1893; Ostrike r 196 4b ; Petroff & Horatschek 2008) 

ri98ll; 



Fischer etalj|2005h . 



1974 


Eriguchi & Susimoto 


iHachisu 


19861: Ansorsetall 



Here we apply the basic idea of the Roche model to rings. 
Thus, we do not choose to have the test fluid rotate in the field of a 
point mass, but in that of a circular line of mass with constant linear 
mass density. In addition to the mass, we thus also have to specify 
the radius of the circle. 

For the comparison of the solutions of the Roche model with 
those of the full problem for polytropes, we use a multi-domain 
spectra l program, much like the one described in lAnsorg et al.l 
(l2003i/) and a similar one tailored to Newtonian bodies with 
toroidal topologies (see lAnsorg & Petroffii2005i for more informa- 
tion). 



2 THE ROCHE MODEL WITH A RING SOURCE 
2.1 Basic Equations 

If we use cylindrical coordinates {g,z,(p), then the gravitational 
potential of a circular line centred on the axis of mass AI and radius 
b reads 



U ■■ 



dip 



GM 

TT y ^Jh^ + q'^ + z"^ - 2bg cos tp 



2GM 



■.K 



(1) 



(2) 



where G is the gravitational constant and K denotes the complete 
elliptic integral of the first kind. 
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Sometimes it will be convenient to use the polar-like coordinates r 
and X defined by 

Q = h — rcosx and z = rsmx. (4) 

For a test fluid rotating uniformly with the angular velocity Q., Eu- 
ler's equation can be written as 



V ( U + h--^'g^ 



0, 



where h is defined by 



h ■- 



p 

f dp' 

J Kp'V 



(5) 



(6) 



and where p is the pressure, /i the mass density and U the potential 
given abovfl We integrate (O and get 



(V) 



* E-mail: S . Horatschek@tpi . uni- jena . de (SH); 
D . Petroff @tpi . uni-jena . de (DP) 



^ For isentropic matter, h is simply the specific enthalpy. 
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where Vb is the constant of integration. Evaluating this equation at 
the surface of the ring, along which the function h vanishes, then 
leads to 



= Vo. 



(8) 



Via this equation, the surface of the ring, described by — rs(x) 
or Zs ~ Zs(f?), and thus the ratio of inner and outer radius and 



A ■- 



Qo 



(9) 



are given implicitly for prescribed M, b, Q and go- Analogously, 
equation l|7j can be used to find surfaces of constant h. 

In what follows, we simplify the equations by introducing the 
dimensionless quantities 



1 

Qo' 



GM' 



(10a) 



(10b) 



r b g z 

r b g z 

U _ Vo _ h _ go 

U ^Vo^h^ GM' 

which implies g\ = A and g^ — 1. In the dimensionless versions 
of the above equations, GM cancels out. Except for two scaling 
constants (e.g. M and £»„), two parameters are necessary to describe 
a ring in the Roche model (e.g. b and f2^). 

2.2 Mass-Shedding Configurations 

Of particular interest are configurations at the mass-shedding limit. 
In this limit, a fluid particle at the outer rim rotates with the Kepler 
frequency, which means that the gravitational force is balanced by 
the centrifugal force alone - the pressure gradient vanishes. For the 
squared angular velocity at the mass-shedding limit we find 



With the potential ^ this gives 



- 



2E{b) 



(11) 



(12) 



where E denotes the complete elliptic integral of the second kind 



(13) 



Contrary to the general case, only one parameter is free, the other 
one is fixed by equation i ll lb . 

Mass-shedding configurations have a cusp at the outer rinfl 
and we will now calculate the associated angle. At the surface, the 
quantity 



V — U 



L ^2 2 

-il g 



(14) 



is constant, see which means that along the surface we have 



= 



dz. 



dg-^ 



= V,,, + 2V,„,-^+V„ 



dg 



dzs 



dg 



^ One can easily show that a configuration with a cusp must be at the mass- 
shedding limit. The converse statement, that mass-shedding configurations 
have a cusp, can be proved if one maJces the very reasonable assumption 
that the pr essure does n ot increase for increasing values of z in the upper 
half-space jPahtj|2007l) . 



where a comma indicates partial differentiation. The equatorial 
symmetry implies that first derivatives with respect to z vanish at 
2 = 0, thus leading to 



g ^ go, z = 



dz.{Q) 
dg 



v.. 



— rj 

^ '-ms ^ ,gg 



(16) 



Hence the full angle is 

5=2 arctan 



02 - TI 



Q—Qo, 2—0 



= 1 



(7,22 

2{l-b'^)E{b) 
{l-P)K{h) ~2E{b)' 



(17) 



cf. Fig. El 



2.3 The Shape of the Surface 

To treat both mass-shedding configurations and the general case, 
we parametrize the angular velocity by 



G [0, 1] 



(18) 



This choice of a arises from the fact that the mass-shedding limit 
(a = 1) poses an upper bound for the angular velocity and a solu- 
tion to l[8) can be found for arbitrarily small Q^. 

Evaluating l[8j at the point {g = go = 1, z = 0) yields 



Vo = 

TT 



aE{b) 



+ 2K{b) 



and at the point {g = gi = A,z = 0) gives 



Vo = 



aA^Ejb) 2^JA 
1-62 5 If, 



Requiring that both expressions for Vb agree leads to 



a{l-A^)E{b) 



2(1-62) 

At an arbitrary surface point {g, Zs,{g)) we get 



(19) 



(20) 



(21) 
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+ zi 



(22) 



Together with J19t this is an implicit equation for the surface func- 
tion Zs{g) or f\{x)- Unfortunately this function cannot be found in 
analytic terms, however it is not difficult to handle it numerically. 
Furthermore, it will be treated analytically for thin rings in the next 
subsection. 

For a prescribed a, the requirement A means that l l21t 
can only be satisfied for 6 larger than some minimal value. It turns 
out that smaller values of 6 no longer describe rings, but spheroidal 
figures. In this case, l l20t and l l21t must be replaced by 

g = 0, z = zj,: Vb = , ^ (23) 



+ z^ 



and 



K{b) 



aE{b) 
2 (1 - 62) ' 



(24) 
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where = Zp/ga denotes the dimensionless polar radius, in other 
words, the ratio of the polar to the equatorial radius. The transition 
from spheroidal to toroidal topologies is described when Zp — oi 
equivalently A = 0. In the limit 6^0, the potential ^ becomes 
that of a point mass A/ at the star's centre and one recovers the 
'standard' Roche model, cf. Appendix lAl 



2.4 Thin Rings 

In the thin ring limit in which A 1, clearly b also tends to I. 
This limit is of particular interest, especially since a self-gravita ting 
thin ring is also tractable to analytical methods ( Kowale wskyl 
1885; Poincare 1885b c d; D vsoiil 1 18921. 1 1 8931 ; [Ostrikej 11964^ : 
Petroff & Horatschelii.200&) . By expanding the surface function 



^-s(x) = E^«(x)(l-6)' 
1=1 

about the thin ring limit 6^1, equation ^ yields 

= q[1 + si(x)cos x] + Insi(x), 

which can be solved by using the Lambert W function, which fulfils 
the equation W{x)e 



(25) 



(26) 



(27) 



a cos X 

For a — > 0, we find si(x) ~ 1 meaning that the cross-section 
tends toward a circle for b ^ 1. For a — 1 (and only for this 
value), si(x) is not differentiable at the point x = ^ although both 
one-sided derivatives exist. One finds 



lim 



dsi(x) 



lim 



dsi(x) 



1, 



(28) 



xTt dx x-lT dx 

which implies that the angle at this point is 5 = n/2, which can 
also be derived from ( 117b in the limit b ^ 1. A series of pictures 
showing the shapes of si(x) for various values of a can be found 
in Fig. [T] Higher terms Si (x) from the series l |25l ) can be found 
iteratively and contain powers of In (l — b). 

To calculate surfaces of constant h, f^(x)i we can easily gen- 
eralize the above calculation. Like {25}, we expand these surfaces 
about the thin ring limit 



(29) 



As mentioned above, the function h vanishes at the surface, i.e. 
fs(x) = ro{x) and Si — ti. Equation (O then gives 



= Q [l + t ? (x) cos x] + In t ^ (x) + tt/i, 



(30) 



and thus 



t?(x)=e- 



-W{ae~ 



cos x)~(^~'^h 



W{ae' 



' cos x) 



a cos X 



(31) 

This function is depicted in Fig.|2]for a = 1 and various values of 
h. For large values of h, we have 



lim ti(x)e 



= 1 



(32) 



meaning that the curves of constant h become circular as the source 
is approached. In the limit a — > we find that 



GM 1 , 

n — ; ill t, 

b TV 



t G (0, 1] (33) 
only depends on the coordinate t :— f/{l — b), but not on x- 




Figure 2. Lines of constant pressure (corresponding to constant h) are 
shown for the mass-shedding ring in the thin ring limit. The function t J (x) 
with Q = 1 is depicted for various values of h in polar coordinates, where 
the axis of rotation is to the left of the cross-section and infinitely far away. 
The surface, which corresponds to = 0, is precisely si(x) and can also 
be found in Fig. [T] The values of h starting from the surface and moving 
inward are 0, 0.001, 0.01, 0.05, 0.1, 0.2, 0.5 and L The value of h tends to 
infinity at the point, where the source is located. 



3 COMPARISON WITH THE FULL PROBLEM 

The comparison of the Roche model with the full problem requires 
the identification of various physical quantities. Whereas AI, Q, 
Vo, A and Zs.{g) are clearly defined in both, the value of b, which 
describes the location of the singular source in the Roche model, 
is not defined in the full problem. There, we choose b to be the 
position of the centre of mass of a meridional cross-section of the 
ring. 

The Roche model contains two scaling parameters, which can 
be 'eliminated' by using the dimensionless quantities introduced in 
l IlOl l, and two physical parameters. On the other hand, for a given 
equation of state, the full problem is determined by specifying a 
scaling parameter and only one physical parameter. When compar- 
ing a specific Roche model with a solution to the full problem, one 
thus has freedom as to how to make such a comparison. One can, 
for example, choose to have b and fl agree and then compare Vo, A 
and the shape of the surface. One could also choose to have b and 
A agree and then compare Q, Vo and the shape of the surface. 

There must be some well-defined way of choosing the addi- 
tional parameter in the Roche model if the full problem is to tend 
to it in some limit. Let us consider rings made up of matter obeying 
the equation of state 

p = K^^'+'^-, (34) 

i.e. polytropes, where K is the polytropic constant and n the poly- 
tropic index. If such rings are expanded about the limit A ^ 1, 
then a solution results, the first few terms of which provide a good 
ap proximation to the full problem over some range of values for 
A tostrike dl 1964bl ; iPetroff & Horatschekll2008h . This range of val- 
ues shrinks to the single point j4 = 1 in the limit n — > oo, i.e. in 
the 'isothermal limit'. This means that the 'isothermal' rings must 
be infinitely thin. Because we expect that only such rings can be 
treated arbitrarily well in the Roche model, we do have a well de- 
fined way of choosing the 'additional parameter': namely A — 1. 
The one remaining free physical parameter, e.g. a, corresponds to 
the single physical parameter one has in the full solution and in- 
terpolates between mass-shed rings and those with circular cross- 
sections. We thus expect that there exist isothermal (and necessarily 
infinitely thin) rings with non-circular cross-sections, cf. Fig.[T] 
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Figure 1. Cross-sections of rings in tlie thin ring limit scaled such that the horizontal extent of each is equal. The function si (x) is depicted for various values 
of a in polar coordinates, where the axis of rotation is to the left of each cross-section and infinitely far away. The dots indicate the origin, where the source of 
the potential is located. 



Those with circular cross-sections were studied in the papers 
mentioned in the last paragraph and the results presented here must 
coincide with those for q = 0. Indeed expression (148 ) for the 

J otentia l and (35) with (151) for the angular velocity in Ostrikeil 
1964bh together with /3"^ > ^ > 1 show that Q^b^ is negligi- 
ble compared to U, as implied by q = and dlSt . Furthermore, 
the logarithmic behaviour J33t can be recovered fro m (120), (122) , 
(125) and (127) in .Ostriker U964bPI (see also iOstrikeilll964j : 
IPetroff & Horatschel 



h,,, = h,-2Kln[ 1 + 



CG [0, 



(35) 



In expression ([33} the function h diverges at the centre t = and 
the surface is located at the finite radius t — 1, whereas in l l35t 
the corresponding function remains finite at the centre ^ = 0, but 
diverges at the surface ^ oo. To compare the expressions for 
h, we thus take a derivative in order to eliminate the physically 
irrelevant constant and consider only the regular portion of the t- 
interval, t > 0, which corresponds to infinite values of ^. This 
leads to 

dh GM 
t— = — (36) 



and 



dt 



^d/iiso 
lim C-JT" 



Trfe 



(37) 



which can indeed be seen to be in agreement upon taking into ac- 
count 

GM = 4-KKb, (38) 

see (146) and the definitions (97) in lPetroff & Horatschekl j2008l) . 
The situation is analogous to that in the spherically symmetric case 
as discussed in Appendix IaI 

If we depart from the thin ring limit, then the comparison with 
numerical solutions to the full problem of a uniformly rotating, 
self-gravitating ring allows us to address the particularly important 
question of how good the Roche model is. Tables[T]and|2]show how 
this model converges toward the full solution for polytropic rings 
as the polytropic index n is increased. The first of these provides a 
comparison of rings with the radius ratio A = 0.7 and the second 
for mass-shedding configurations (q = 1 for the Roche model). In 
both cases, the value of b for the Roche model was chosen to agree 
with the numerical one. The second column, k, gives an indication 
of how concentrated the mass is. If we consider the cross-section 
of a ring, and define the two points in the equatorial plane at which 
the density falls off to half of its maximal value 



^(£la,2 = 0) = ^J,igb,z = 0) 



Qb > £>a 



(39) 



Table 1. Comparison between polytropic rings with A = 0.7 and the cor- 
responding Roche configuration with the same 5. A denotes the difference 
between the values of the Roche configuration and the values of the full 
solution. 



n 


ft 




b 




/l^nuni 


Ayo/Vonum 


1.0 


6.3 X 10- 


1 


8.48 X 10- 


1 


2.6 X IQ- 


.2 


7.4 X 10-3 


1.5 


4.9 X 10" 


1 


8.46 X 10" 


1 


1.9 X IQ- 


.2 


6.0 X 10-3 


3.0 


2.7 X 10" 


1 


8.42 X 10- 


1 


8.2 X IQ- 


'3 


3.1 X 10-3 


5.0 


1.3 X 10- 


1 


8.34 X 10- 


1 


2.6 X IQ- 


'3 


1.2 X 10-3 


7.0 


7.2 X 10" 


2 


8.26 X 10- 


1 


8.6 X 10- 


-4 


4.7 X 10-* 



then K is defined to be the ratio of the distance between these two 
points to the total width of the ring's cross-section 



(40) 



3 Equation (122) in lOstrikeJ h964d) should read p = Ae-* 



Only for K ^ 1 is the Roche model expected to give good results. 
We find that k < 0.2 for polytropic rings with n > 3. Such poly- 
tropes always have values for the radius ratio A > 0.45, which 
implies that Roche models with a smaller radius ratio do not pro- 
vide a particularly close approximation to any toroidal polytrope. 

The shapes of cross-sections of mass-shedding rings for var- 
ious polytropic indices are plotted in Fig. |3] Here, the differences 
that arise due to different prescriptions of the parameters are evi- 
dent for small values of the polytropic index. When A is chosen to 
coincide between the Roche model and the full solution, then the 
two surfaces do not differ appreciably, even for n — 1. The ra- 
dius ratio itself, and consequently the shape, is quite different for a 
Roche model with b as prescribed from the solution to the full prob- 
lem with n — 1. These differences vanish as n is increased and the 
two prescriptive choices converge to the full solution as they must. 

Fig. |4] provides an interesting comparison between mass- 
shedding configurations for polytropes and the Roche model. The 
angle S of the cusp as given by l ll7t is plotted over the whole in- 
terval b £ [0, 1] as a solid line. This curve is compared to two 
distinct polytropic mass-shedding sequences generated by varying 
the polytropic index n. The dotted line describes rings and merges 
together with that of the Roche model in the thin ring limit, i.e. for 
6-^1. One the other hand, the dashed line describes (spheroidal) 
stars and merges together with that of the Roche model for 6^0. 
The change in topologies in the Roche model from toroidal to 
spheroidal takes place for fe ~ 0.56 as can be calculated by solving 
for b in l|2T) with ^ = and a = 1. 

In summary, we can say that as with stars, the Roche model 
presented here for rings is seen to yield a very good approximation 
to the full problem in many cases. Moreover, it presumably pro- 
vides exact results in the appropriate limit and offers new solutions 
that had not been found using perturbative approaches. 
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Table 2. Comparison between polytropic rings in mass-shedding and the corresponding Roche configiu-ation with the same b. A denotes the difference between 
the values of the Roche configuration and the values of the fuU solution. 



n K b AA/Annm ACl/Clnum AS/Snam AVo/Vq 



1.0 


5.3 X IQ- 


1 


5.77 X 10" 


1 


-3.7 X 10- 


-1 


-3.8 X 10" 


-2 


2.7 X 10" 


-2 


-4.1 X 10- 


-2 


1.5 


4.0 X 10" 


1 


6.04 X 10" 


1 


-1.0 X 10" 


-1 


-2.0 X 10" 


-2 


1.3 X 10" 


-2 


-2.1 X 10- 


-2 


3.0 


2.1 X 10" 


-1 


6.53 X 10" 


-1 


-1.2 X 10" 


-2 


-3.9 X 10" 


-3 


2.1 X 10" 


-3 


-4.4 X 10" 


-3 


5.0 


1.1 X 10" 


-1 


6.92 X 10" 


-1 


-2.0 X 10" 


-3 


-7.7 X 10" 


-4 


3.4 X 10" 


-4 


-8.8 X 10" 


-4 


7.0 


5.7 X 10" 


-2 


7.20 X 10" 


-1 


-4.7 X 10" 


-4 


-2.0 X 10" 


-4 


7.8 X 10" 


-6 


-2.4 X 10" 


-4 
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3.0 



5.0 



7.0 
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Figure 3. Comparison between the cross-sections of polytropic rings at the 
mass-shedding limit (dotted Unes) and those of the corresponding Roche 
configurations, i.e. with a = 1 (solid lines). In the left panels the Roche 
model was chosen such that b agrees with the numerical value, whereas in 
the right panels, the values of A were chosen to coincide. The ticks in the 
'centre' indicate the values of b. When A is prescribed (right panels), it 
turns out that 6„um < feRoche- 



120° 



110° 



100° - 



90 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 4. The angle 5 of the cusp for the mass-shedding configurations in 
the Roche model (solid line) is compared with the numerical result for poly- 
tropic rings (dotted line) and polytropic spheroidal configurations (dashed 
line). The larger dots denote the ring configurations with n = 1/2, 1, 3/2, 
3, 5 and n = 7, where smaller values of n correspond to smaller values of 6, 
and the white squares denote the spheroidal configurations with n = 1/2, 
1, 2, 3 and n = 4, where larger values of n correspond to smaller values 
of 5. In the thin ring limit, where we expect that the Roche model provides 
arbitrarily good results to the full problem for 'isothermal' rings (n — > oo), 
wefind<5 -> 90°. 
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APPENDIX A: THE STANDARD ROCHE MODEL 

As mentioned in the intro duction, the R o che model with a 1/r- 

J otent i al was discussed by 
197lh : IShapiro & Shibatal 



where 



Rochd ( Il873h : IZel'dovich & Novikovl 
2002*) and provides a unique surface 
function describing a mass-sh edding star. To derive this function, 
we follow lMeinel et alj j2008l) and begin by writing down equation 
([8} with the potential Us of a point mass 



Vo 



GM 



-il g — 0. 

2 ^ 



(Al) 



By considering the point g — 0, we can relate the constant Vo to 
the polar radius Zp 

GM 



Q—Qi.1 ,2 — 



z 

For mass-shedding stars, the relation 

dU_ 

dg 

where g^, is the equatorial radius, leads to 
GM 

and the surface equation JAU becomes 



go^l 



g„il , 



= 1. 



(A2) 



(A3) 



(A4) 



(A5) 



Evaluating this expression at the equator, 2 = 0, tells us that for 
mass-shedding fluids in the Roche model, the radius ratio 

z„ 2 



£-0 



(A6) 



follows. With this relationship, the curve describing the fluid's sur- 
face can be rewritten as 



It then follows that 



lim ^ = -^3, 
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which means the interior angle of the mass-shedding cusp is 120° 
(see Fig.Hfor h 0). 

In general, the Roche model is expected to approach the full 
solution when the matter is arbitrarily concentrated. For static poly- 
tropes, the full problem leads to the Lane-Emden equation, which 
turns out to describe arbitrarily concentrated matter if n = 5. We 
shall now show the agreement between these stars and the results 
given by the Roche model. 

In the static case, the Roche model gives 



h = GM \ ^ - 4- 
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and thus 



The solution to the full problem is 

^ ^2\-5/2 

MS = ' ' 

and thus 



e 



?e [o,oo) 



(A9) 



(AlO) 



(All) 



I' 



3K 



2TTGfJ., 



4/5 



(A13) 



is the dimensionless radius and x the true radius. The surface ex- 
tends out to infinity and since the regular portion of the 7?-interval, 
R > 0, corresponds to infinite values of the radius x, we consider 

2d/l5 



lim x^^^^ = —6\/3Klij,l^^ . 
a;— >oo da; 
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Calculating the mass for the full problem, 

oo 

GM = 4ttG J fisx'^ dx = 47r^/3^Vc = eVSKlfij^^, (A15) 



shows the equivalence of dAlOt and ( IA14t . 
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